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Scaled Brownian motion (SBM) is widely used to model anomalous diffusion of passive tracers 
in complex and biological systems. It is a highly non-stationary process governed by the Langevin 
equation for Brownian motion, however, with a power-law time dependence of the noise strength. 
Here we study the aging properties of SBM for both unconfined and confined motion. Specifically, 
we derive the ensemble and time averaged mean squared displacements and analyze their behavior 
in the regimes of weak, intermediate, and strong aging. A very rich behavior is revealed for confined 
aging SBM depending on different aging times and whether the process is sub- or superdiffusive. We 
demonstrate that the information on the aging factorizes with respect to the lag time and exhibits 
a functional form, that is identical to the aging behavior of scale free continuous time random walk 
processes. While SBM exhibits a disparity between ensemble and time averaged observables and is 
thus weakly non-ergodic, strong aging is shown to effect a convergence of the ensemble and time 
averaged mean squared displacement. Finally, we derive the density of first passage times in the 
semi-infinite domain that features a crossover defined by the aging time. 

PACS numbers: 05.40.-a 


I. INTRODUCTION 

Deviations from normal Brownian motion were re¬ 
ported already in the work of Richardson on the spread¬ 
ing of tracers in turbulent flows [lj, and deviations from 
the Brownian law are discussed by Freundlich and Kruger 
[0]. Today anomalous diffusion is typically defined in 
terms of the power-law form 

(1) 

of the mean squared displacement (MSD) 0 0. De¬ 
pending on the value of the anomalous diffusion exponent 
we distinguish the regimes of subdiffusion (0 < a < 1) 
and superdiffusion (a > 1), including the special cases of 
Brownian motion ( a = 1) and ballistic transport (a = 2). 
The generalized diffusion coefficient K* a in Eq. © has the 
physical dimension cm 2 /sec a . 

Anomalous diffusion is observed in a wide range of sys¬ 
tems, including fields as diverse as charge carrier motion 
in amorphous and polymeric semiconductors Jjl 0, dis¬ 
persion of chemicals in groundwater aquifers [7[, particle 
dispersion in colloidal glasses (§}, or the motion of tracers 
in weakly chaotic systems 0. With the rise of experi¬ 
mental techniques such as fluorescence correlation spec¬ 
troscopy or advanced single particle tracking methods, 
the discovery of anomalous diffusion has gone through 
a sharp rise for the motion of endo gen ous and artificial 
tracers in living biological cells flOl 1151 . Concurrent to 
this development an increasing amount of anomalous dif¬ 
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fusion studies is reported in artificially crowded environ¬ 
ments mimicking aspects of the superdense state of the 
cellular fluid [13, [l8| . Within and along lipid membranes 
anomalous diffusion was found from experiment and sim¬ 
ulations El- 

Brownian motion is intimately connected with the 
Gaussian probability density function describing the spa¬ 
tial spreading of a test particle as function of time. This 
Gaussian is effected a forteriori by the central limit the¬ 
orem, as Brownian motion is well described on a stochas¬ 
tic level by the Wiener process. Anomalous diffusion 
loses this universal character, and instead different sce¬ 
narios corresponding to the physical setting need to be 
considered. Among the most popular models we men¬ 
tion the Scher-Montroll continuous time random walk 
(CTRW), in which individual jumps are separated by 
independent, random waiting times 0E3. If the dis¬ 
tribution of these waiting times is scale free, subdiffusion 
emerges [13 ■ Fractional Brownian motion and the closely 
related fractional Langevin equation motion are stochas¬ 
tic processes fueled by Gaussian yet power-law correlated 
noise [2ll , [22| . Anomalous diffusion emerges when a con¬ 
ventional random walker is confined to move on a ma¬ 
trix with a fractal dimension 0 E Hi- Stochastic 
processes with multiplicative noise, corresponding to a 
space-dependent diffusion coefficient, also effect anoma¬ 
lous diffusion [25i . [26| . A contemporary summary of dif¬ 
ferent anomalous diffusion processes exceeding the scope 
of this introduction is provided in Ref. E ■ 

Here we deal with the remaining of these popular 
anomalous diffusion models, namely Scaled Brownian 
Motion (SBM). SBM is a highly non-stationary process 
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defined in terms of the stochastic equation 

^ x m, (2) 

which is driven by white Gaussian noise of zero mean 
(£(i)) = 0 and with autocorrelation (£(ti)£(t 2 )} = d(t\ — 
£ 2 ) • The explicitly time dependent diffusion coefficient 
is taken as 

J(r{t) = aK* a t a ~\ (3) 

We allow a to range in the interval (0, 2), such that the 
process describes both subdiffusion and sub-ballistic su¬ 
perdiffusion. The idea of a power-law time dependent 
diffusion coefficient is essentially dating back to Batche¬ 
lor (albeit he used a = 3) in his approach to Richardson 
turbulent diffusion [28], SBM, especially in its subdiffu- 
sive form, is widely used to describe anomalous diffusion 
[29{ . SBM was studied systematically in Refs. [26ll30l432j | . 

In stationary systems correlations measured between 
two times t\ and £2 are typically solely functions of the 
time difference, f(\t± — £ 2 !)- In non-stationary systems 
this functional dependence is generally more involved, 
e.g., it can acquire the form ffo/ti) [231 ■ 1 11 such a non * 
stationary setting the origin of time can no longer be 
chosen arbitrarily. This raises the question of aging, that 
is, the explicit dependence of physical observables on the 
time span t a between the original preparation of the sys¬ 
tem and the start of the recording of data. Traditionally, 
aging is considered a key property of glassy systems [3J] . 
The aging time t a can be adjusted deliberately in certain 
experiments, such as for the time of flight measurements 
of charge carriers in polymeric semiconductors in which 
the system is prepared by knocking out the charge carri¬ 
ers by a light flash [gj. Similarly, aging could be checked 
directly in blinking quantum dot systems, in which the 
initiation time is given by the first exposure of the quan¬ 
tum dot to the laser light source. In other systems, for 
instance, the motion of tracers in living biological cells, 
the aging time is not always precisely defined. In such 
cases it is therefore important to have cognisance of the 
functional effects of aging as developed here. 

In the following, we will analyze in detail the aging 
properties encoded in the SBM dynamics in both un¬ 
confined and confined settings. For free aging SBM in 
Section [TT] we show that the result for the time averaged 
MSD factorizes into a term containing all the information 
on the aging time t a and another capturing the physi¬ 
cally relevant dependence on the lag time A. This fac¬ 
torization is identical to that of heterogeneous diffusion 
processes and scale-free, subdiffusive CTRW processes. 
In Section Mil we explore the aging dynamics of confined 
SBM. For increasing aging time t a it is demonstrated that 
the non-stationary features of SBM under confinement 
are progressively washed out, a feature, which is impor¬ 
tant for the evaluation of measured time series. Section 
llVI reuorts the first passage time density on a semi-infinite 
domain for aged SBM which includes a crossover between 
two scaling regimes as a result of the additional time scale 
introduced by t a ■ Finally, Sec. f\~l concludes this paper. 


II. AGEING EFFECT ON UNCONFINED SBM 


The position autocorrelation function (covaria nce) for 
SBM in the conventional (ensemble) sense reads [30j 

(x{ti)x(t 2 )) = 2AT* min(ti, t 2 )“- (4) 

For an aged process, in which we measure the MSD start¬ 
ing from the aging time t a until time t, the result for the 
MSD thus becomes 


{x 2 {t))a = (\x(t a +t) -x(t a )] 2 ) 

= 2 K*[(t+t a r~a. ( 5 ) 


For a non-aged process with t a = 0 the standard scaling 
© of the MSD is recovered, as it should. In the aged 
process, the MSD © is reduced by the amount accu¬ 
mulated until time t a , at which the measurement starts. 
The limiting cases of expression © interestingly reveal 
the crossover behavior 


(x 2 (t)) 


2aK*t%~ 1 t, t„»i 
2 K* a t a , t > t a 


( 6 ) 


While for weak aging (t a <C t) the aged MSD © be¬ 
comes identical to the non-aged form ©, for strong ag¬ 
ing (t a t) the scaling with the process time t is linear 
and thus, deceivingly, identical to that of normal Brow¬ 
nian diffusion. However, the presence of the power t 
is reminiscent of the anomaly a of the process. We note 
that the behavior © and thus © is identical to the re¬ 
sult for the subdiffusive CTRW [35l. [36} as well as aged 
heterogeneous diffusion processes with a power-law form 
of the position dependent diffusivity [37j . 

In single particle tracking experiments [38l - [4l| one 
measures the time series x(t) of the position of a labeled 
particle, which is then typically evaluated in terms of the 
time averaged MSD. For an aged process originally ini¬ 
tiated at t = 0 and measured from t a for the duration 
(measurement time) t this time averaged MSD is defined 
in the form [?!(• 


*2(A) = 


t - A 


r 


t + ta-A r 2 

x(t' + A) — x(t') dt', (7) 


as a function of the lag time A and the aging time t a . 
Averaging over an ensemble of N individual trajectories 
in the form 



1 N _ 


( 8 ) 


the structure function ([x(t' + A) — x(t')} 2 ) in the inte¬ 
gral of expression © can be evaluated in terms of the 
covariance ©• The exact result reads 



2 K* 

(a + l)(i-A) 


(t + t a ) a+1 


~(t + t a — A)“ +1 + f“ +1 


(t a + A)“ +1 


( 9 ) 
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In the absence of aging, we recover the known result [26j, 

mm 

(j*(A))~2 K*-^. (10) 

Its linear lag time dependence contrasts the power-law 
form of the ensemble averaged MSD m and thus demon¬ 
strates that this process is weakly non-ergodic in the 
sense of the disparity [13, IH, EH 

(nA)}^<* 2 (A)>. (11) 

The equivalence and therefore ergodicity in the Boltz¬ 
mann sense is only restored in the Brownian case a = 1. 
In the presence of aging, expansion of expression © in 
the limit t,t a ^> A of short lag times yields 

(sm) - x (P{A)) , (12) 

in which we defined the so-called aging depression as 

A a (z) = (1 + z) a - (z) a . (13) 

In this experimentally relevant limit all the information 
on the age of the process is thus contained in the aging 
depression A a , and the physically important dependence 
on the lag time A factorizes, such that Eq. (fl^l) con¬ 
tains the non-aged form m- Result m is identical to 
the behavior of aged subdiffusive CTRW [36j and het¬ 
erogeneous diffusion processes In the limit t a t 
of strong aging, the time averaged MSD © remarkably 
reduces to the form 

(WAj) = 2 (14) 

In this limit, the time averaged MSD thus becomes equiv¬ 
alent to the aged ensemble averaged MSD, (d 2 (A)) = 
(x 2 (A)) a , as evidenced by comparison with result ([5]). 
In this limit, that is, ergodicity is restored, as already 
observed for aged CTRW processes [3Cj ] . 

Figure [T] shows the behavior of the ensemble and time 
averaged MSD for unconfined SBM at different aging 
times in the subdiffusive case with a = 1/2. The thin 
lines depict the simulations results for the time averaged 
MSD for 20 individual trajectories. The first observa¬ 
tion is that the amplitude spread between these 20 time 
traces is fairly small. Note that the larger scatter for 
longer lag times A is due to worsening statistics when A 
approaches the trace length t. The circles in Fig.[T]corre- 
spond to the average over the 20 different results for the 
time averaged MSD. The latter compare very nicely with 
the theoretical expectation m ■ Finally, the thick green 
line is the theoretical result ([5]) for the ensemble averaged 
MSD. The detailed behavior in the three different aging 
regimes is as follows: 

(i) In the non-aged case ( t a = 0, top panel of Fig. [TJ 
the power-law growth (x 2 (t)) ~ t a of the MSD contrasts 
the linear form (<5 2 (A)) ~ A, this disparity being at the 
heart of the weak ergodicity breaking [2a, E3, |33 ■ 



t, A 



FIG. 1: Ensemble and time averaged MSD for SBM with 
a = 1/2. Thin lines: time averaged MSD for 20 individual 
trajectories from simulations of the SBM Langevin equation 
m with trajectory length t = 10 5 . Circles: averages over 
those 20 trajectories. Black thin line: theory result m- 
Thick green line: ensemble averaged MSD Q. Three dif¬ 
ferent aging times were considered (top to bottom): (a) non¬ 
aged case t a = 0, (b) weak aging case t a = 10 3 , and (c) strong 
aging case t a = 10 6 . In all simulations A/ = 1/2. 


(ii) In the weak aging case ( t a = 10 2 , middle panel 
of Fig. [TJ a major change is visible in the behavior of 
the MSD, namely, we see the crossover from the aging- 
dominated linear scaling (x 2 (t)) ~ to the anoma¬ 

lous scaling (x 2 (t)) ~ t a , encoded in Eqs. ©. The be¬ 
havior of the time averaged MSD is largely unchanged in 































4 



FIG. 2: Ensemble and time averaged MSD for SBM ([2j with 
a = 3/2 in the strong aging case, t a = 10 6 . The observation 
time is t = 10 5 . The spread of the 20 single trajectory time 
averages is fairly small. As before, the ensemble and time 
averaged MSDs coincide, apparently restoring ergodicity. 


comparison to case (i). 

(iii) In the strong aging case (t a = 10 6 , bottom panel 
of Fig. [lj we see the apparent restoration of ergodicity: 
ensemble and time averaged MSDs coincide, as given by 
Eq. (TUl) . 

The convergence of the ensemble and time averaged 
MSDs in the strong aging case for the superdiffusive case 
with a = 3/2 is nicely corroborated in Fig. [2] 


III. AGEING EFFECT ON CONFINED SBM 

In many cases an observed particle cannot be consid¬ 
ered free during the observation. Examples contain par¬ 
ticles moving in confined space, for instance, within the 
confines of living biological cells [3, 0]. Similarly, par¬ 
ticles measured in optical tweezers setups experience a 
confining Hookean force [3 [3 EH- As a generic ex¬ 
ample for confined SBM we consider the linear restoring 
force — kx(t ) with force constant k. The corresponding 
stochastic equation for this confined SBM reads (32| 

= -kx(t) + ^2aK **(“-!) x £(f), (15) 

where, as before, £(f) represents white Gaussian noise of 
zero mean. The covariance in this confined case yields in 
the form [HI 

(s(ii)a;(i 2 )) = 2 7F*f“ e - fc(tl+t2) M(a, a + 1,2 kt x ) (16) 

for ti < t 2 in terms of the confluent hypergeometric func¬ 
tion of the first kind, also referred to as the Kummer 
function (32, 47j. Based on this result we now present 
the ensemble and time averaged MSDs. 


A. Ensemble averaged MSD of confined SBM 

The ensemble averaged MSD for aging SBM, 
{x 2 {t)) a = ([x(t a +t) - x{t a )] 2 } becomes 

(x 2 (t)) a = 2Jt x {t a +t) + 2y/Ki(t a ) - 4e- fet .#i (t a ), (17) 

where we used the abbreviation 

(t) = K*t a exp(—2 kt)M(a, a + 1,2 kt). (18) 

In the limit k —> 0 of vanishing confinement, Eq. © for 
free SBM is readily recovered from the property M(a, a+ 
1 , 0 ) = 1 . 

We now discuss the result (ED in the three limits of 
the non-aged, weakly aged, and strongly aged processes. 
The analysis reveals a rich behavior depending on the 
values of the aging time t a and the anomalous diffusion 
exponent a. For sub- and superdiffusion, respectively, 
the various crossovers are displayed in Figs. [3] and [4] 

(i) In the absence of aging (t a = 0), we get back to the 
result 

(: x 2 {t)) = 2Jt x (t) (19) 

reported in Ref. [3]. For t ^ 1/k this reduces to the 
non-aged free SBM result Q, while in the long time limit 
t 1/k we use the expansion 

M(a,a + l,z)~a^^ (20) 

z 

of the Kummer function to obtain [32f 

nK* 

(x(t) 2 ) ~ (21) 

This result underlines the inherently non-stationary char¬ 
acter of SBM: for subdiffusion the MSD ( x(t ) 2 ) progres¬ 
sively decays, while for superdiffusion it increases. This 
property reflects the time dependence of the temperature 
encoded in the diffusivity © j32| • This non-aged behav¬ 
ior is shown in Figs. [3] and [4] as the grey lines for two 
different strengths k of the external confining potential. 
How does aging modify this behavior? 

(ii) We first consider the case C 1/fc. If in addition 

t <C 1/k, this is but the above result © for free aging 
SBM. However, care needs to be taken when t 1/k. 

From Eq. ED, we find that the first two terms (the third 
one is exponentially small in t and can be neglected) lead 
to the asymptotic behavior 

rvK* 

(x 2 (t)) a ~ n ^t‘*- 1 +2K*X. (22) 

This implies that for subdiffusion (0 < a < 1) the first 
term tends to zero and the leading behavior is the plateau 

{x 2 {i)) a ~2K* a r a . (23) 

Even for very weak aging, the ensemble averaged MSD 
( x 2 (t)) a becomes independent. When experimental data 
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FIG. 3: Ensemble averaged MSD ( x 2 (t )) for confined aging SBM in the subdiffusive case with a = 0.5 at different aging times: 
(i) non-aged ( t a = 0) denoted by the grey lines; (ii) weakly aged ( t a = 0.1) denoted by the blue lines; and (iii) strongly aged 
( t a = 10 6 ) denoted by the orange lines. In all cases, the full lines correspond to the force constant k = 0.1, while the dashed 
lines stand for k = 0.01. The green line at the bottom of the graph is a blowup (( x 2 (t )) 4 of the case t a = 10 6 and k = 0.1) in 
which the crossover between the two plateaux is better visible. 


are evaluated and the exact equivalence t a = 0 is not 
guaranteed, the erroneous conclusion could be drawn 
that the process is stationary. Note, however, that result 
(l23l) is independent of the strength k of the confining po¬ 
tential and only depends on the diffusion coefficient IT* 
and the aging time t a , mirroring the fact that this term 
stems from the initial free motion during the aging pe¬ 
riod. Conversely, for superdiffusion (a > 1) the leading 
order term indeed shows the growth 

rvK* 

(x 2 (t)) a ~ (24) 

k 

of the ensemble averaged MSD. The weakly aged behav¬ 
ior is shown in Figs. [3] and H] as the blue lines. 

(iii) With the asymptotic expansion (1201) of the Rum¬ 
mer function, we find that in the strong aging regime 
t a 1/k the ensemble averaged MSD becomes 

{x 2 {t))a ~ ak~ x K* a [(t a + t) a ~ x + (1 - 2e" fet )] . 

(25) 

At short times t <C 1/A;, this leads us back to the uncon¬ 
fined result ( x 2 (t)) a ~ 2a/\*f“ _1 f of Eq. (0. At long 
time t 1/A;, however, we have to distinguish two dif¬ 
ferent regimes. First, for t a 3> t 1/k we obtain the 
plateau 

(x 2 (t))a ~ (26) 

which differs from the above result (E3D by the factor 
of two. Second, for t S> t a 1/A; the leading order 


according to Eq. C5l) again differs between sub- and su- 
perdiffusive motion. For 0 < a < 1 the plateau 

(x 2 (t)) a ~ (27) 

emerges. Note, however, that in comparison to Eq. (1261) 
we now have half the amplitude. In the superdiffusive 
case a > 1 we recover result CHD. This intricate behavior 
is shown in Figs. [3] and []] as the orange lines. In Fig. [3] 
we pronounce the crossover between the two plateaux by 
plotting the fourth power of the ensemble MSD as the 
green line. 


B. Time averaged MSD of confined SBM 


The time averaged MSD for confined SBM can be de¬ 
rived by substituting the above covariance m into the 
integral 0. By help of the relation [46| 

J y a e~ v M(a, 1 + a, y)dy = 

— x 1+a e- x M(l+a,2 + a,x) (28) 
1 + a 

this procedure yields the general result 
2 K* 


*2(A)) = 


(t — A)(l + a) - 

(t + t a — A) — ^2(ta) 


+ ta) ~ (ta + A) 


—2e 


-k A 


^.#2(A + t a — A) — ^#2(Aa)^ ,(29) 
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FIG. 4: Ensemble averaged MSD ( x 2 (t )) for confined aging SBM in the superdiffusive case with a = 1.5 at different aging 
times: (i) non-aged (t a = 0) denoted by the grey lines; (ii) weakly aged ( t a = 0.1) denoted by the blue lines; and (iii) strongly 
aged (t a = 10 6 ) denoted by the orange lines. In all cases, the full lines correspond to the force constant k = 0.1, while the 
dashed lines stand for k = 0.01. In all cases the terminal scaling ~ t a ~ l is reached. 


where we used the abbreviation 

= t 1+a e~ 2kt M (1 + a,2 + a, 2kt). (30) 

(i) In the limit fc —I 0 we recover the result <[9l) of 
unconfined aging SBM, while the complete absence of 
aging restores the result from Ref. [H]. In the presence 
of confinement, we distinguish the following regimes. 

(ii) We now consider the case when the aging time 
is short compared to the relaxation time of the system, 
t a 1 /k. From the general expression (T2IJ1) we then find 
the following behaviors: (a) when in addition the lag time 
is short (f > 1/1: > A > f a ) we recover the non-aged 
result (1101) with its linear scaling in the lag time A. (b) 
When the lag time is long (f> A>l/fc> t a ), we find 
the plateau 


(A Q 2 (A)) ~ (31) 

known from the non-aged case [27[. (c) Finally, when the 
lag time approaches the length t of the time series, the 
time averaged MSD 

(W*)) ~ (32) 

becomes equivalent to the ensemble averaged MSD, 

Eq. (l24j) . In contrast to the ensemble averaged MSD, 

we thus find that the time averaged MSD is not affected 
by short aging times as compared to the relaxation time 
scale, t a <C 1 /k. 


(iii) The second, more interesting case corresponds to 
long aging times compared to the relaxation time scale, 
t a 1/k. When also t 1/k, the result is indepen¬ 
dent of the specific magnitude of the lag time. From the 
general expression (E^l) by help of relation EU1) we obtain 



K* 


[t + t a ) a — (A + t a ) C 


k(t — A) 

+(1 — 2e -fcA )[(i + t a ~ A)“ — t“]l .(33) 


If we now consider the regime in which the lag time is 
short, t, t a ^ 1/k A, we obtain result (fl2l) with the 
aging depression (fl3l) from unconfined aging SBM. In the 
opposite limit t,t a ^> A 1/k when the lag time is long 
compared to the relaxation time, we find 

(%(£))-Aa{ta/t)(s*(*)), (34) 

where ^<5 2 (A)^> is equal to expression (1311) and A a (z) is 

again the aging depression (fl3l) . In the strong aging limit 
t, t a 3> A » 1/k, that is, the aged time averaged MSD is 

generally given by (^(A)^) - A a (t a /t) ^8 2 (A)^ for any 

lag time. Similar to subdiffusive CTRW processes J2fj, 
the occurrence of the factor A a appears like a general 
feature for the aging dynamics of SBM. 

Figure [5] shows the behavior of the ensemble and time 
averaged MSD for confined SBM at different degrees of 
aging. The graphs represent the full behavior according 
to Eqs. (fT71) and (EH1) . In the absence of aging, the ini¬ 
tial linear growth (S 2 ( A)) ~ A of the time averaged MSD 
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FIG. 5: Ensemble and time averaged MSD for confined SBM 
for a = 1/2. From top to bottom, the panels represent the 
non-aged (t a = 0) case, the case of weak aging ( t a = ICG 1 ), 
and the case of strong aging ( t a = 10 6 ), where the observation 
time is chosen as t = 5 x 10 4 . The lines represent Eqs. G3 
and (HU). The force constants k are indicated in the panels. 
Note that the time averaged MSD indeed converges to the 
ensemble MSD in the limit A —¥ t, compare the discussion in 

Ref. Un¬ 


crosses over to an apparent plateau, contrasting the func¬ 
tional behavior of the ensemble average: at short times, 
we observe the power-law growth ( x 2 (t)) a ~ t a of un¬ 
confined SBM, while after engaging with the confining 
potential, the monotonic decrease (x 2 (t)) a — t a ~ l re¬ 
flects the temporal decay of the noise strength (i.e., the 


xio 3 



FIG. 6: Ensemble and time averaged MSD for confined SBM 
with a = 3/2 in the strong aging case with t a = 10 6 . The 
observation time is t = 5 x 10 4 . 

temperature) (32j. When aging effects come into play, 
the ensemble averaged MSD displays notable differences. 
In the case of weak aging displayed in the middle panel 
of Fig. [5] deviations from the power-law decay of (x 2 (t)) a 
become apparent for longer times, t t a 1/k. Even¬ 
tually the convergence to a common value independent of 
the force constants is observed, as predicted by Eq. ([251) . 
Finally, in the strong aging limit, the ensemble and time 
averaged MSD ar e equivalent and ergodicity is seemingly 
restored: (<5^(A)) = ( x 2 (A)) a , as can be witnessed in the 
bottom panel of Fig. [5j The apparent equivalence of en¬ 
semble and time averaged MSDs in the strong aging limit 
is also proven in the superdiffusive case for a = 3/2 in 
Fig. [6] The slight discrepancy remaining between time 
and ensemble averaged MSD in the latter strong (but fi¬ 
nite) aging case is shown in Fig. (0. This Figure also 
demonstrates the convergence of ensemble and time av¬ 
eraged MSDs. 

IV. FIRST PASSAGE TIME DENSITY 

Apart from the MSD the first passage behavior is a sig¬ 
nature quantity of a stochastic process. We here study 
how aging changes the first passage statistic of SBM in 
the semi-infinite domain. The probability density func¬ 
tion (PDF) of first passage is found by solving the SBM 
diffusion equation with the time dependent coefficient 
JT(t), 

d o2 

—P(x,t)=Jfr(t)—P(x,t), (35) 

however, with the aged initial condition 

^ (36) 

This aged initial condition emerges from a 6(x) peak for 
a system initialized some aging time t a before. In this 
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FIG. 7: Full behavior of aging confined SBM for t a = 10 6 with a = 1/2 (left) and a = 3/2 (right), demonstrating the 
convergence of the time averaged MSD to the ensemble averaged MSD in the limit A —> t, shown for two different potential 
strengths, as indicated in the panels. 


setup t measures the time span from the aged initial con¬ 
dition (13611 . To obtain the first passage PDF for the semi¬ 
infinite domain we solve the SBM diffusion equation (1351) 
for unconfined motion with the aged initial condition (1361) 
and then use the method of images. For the PDF of the 
aged process we obtain 


P{x,t) 


___ e ~ x2 / 4K a (*“+*“) 

y/4TvK*(t% +t a ) 


(37) 


In the presence of an absorbing boundary at the origin, 
the survival probability for a process initiated originally 
in xq > 0 is therefore given by 


/ [P(x — xo, t) — P(x + xq, f)] dx. (38) 

Jo 

Substituting the aged PDF (l3f>l) into this expression 
yields 


For a = 1 (Brownian motion) and in the absence of aging 
( t a = 0) we recover the well known Levy-Smirnov distri¬ 
bution. Result pup exhibits a crossover relative to the 
aging time, 


t a 3a/2 t a 1 , t a > t, (Xq/P-K*]) 1 /" 

t-'~ a/2 , i»t«,(*o/[4tf£]) 1/a 

(41) 

In the strong aging limit the scaling exponent is —(1 — a), 
and we observe the explicit presence of the aging time t a 
with exponent a. For weak aging, the scaling exponent of 
t is —(1 + a/2), as known from subdiffusive CTRW pro¬ 
cesses. However, the detailed crossover behavior is differ¬ 
ent, compare Ref. j48j. We also note that for fractional 
Brownian motion the anomalous diffusion exponent a en¬ 
ters oppositely [30j, [Hj]. Fig. © shows the crossover of 
the first passage density for aging SBM. 


P - 


axo 




y{t) = erf 


x 0 


\/4 


(39) 


V. CONCLUSIONS 


in terms of the error function. The first passage PDF 
follows from the relation p(t) = —dy{t)/dt, 


pit) 


axo t° 


\J 4nK * (t“ + t a )' 


: exp 


x 0 \ 

4 K*(t2 + t«)J ' 

(40) 


SBM is possibly the simplest anomalous diffusion 
model, and it is therefore widely used in literature. De¬ 
spite its apparent simplicity SBM exhibits weak ergodic- 
ity breaking in the sense that we observe a distinct dis¬ 
parity between the ensemble and time averaged MSDs 
of this process [2(1 [53, [31/ |32|. It is therefore a natu¬ 
ral question to explore the aging effects of SBM, i.e., the 
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FIG. 8: First passage time density p(t) for a = 1/2, with 
*o = l, and aging time t a = 10 6 . As shown by the dashed 
line, the crossover between the aging dominated slope —1/2 
to the slope —5/4 is distinct. 

explicit dependence of physical observables on the time 
span t a between the original system preparation. We 
here showed how the ensemble and time averaged MSDs 
depend on t a for both the unconfined and confined cases. 

For unconfined aging SBM we obtained the exact de¬ 
pendence on the aging time t a and observed a strik¬ 
ing similarity to both subdiffusive CTRW with scale-free 
waiting time distributions and heterogeneous diffusion 
processes with power-law position dependence of the dif- 
fusivity. In particular, for short lag times the time av¬ 
eraged MSD factorizes into the non-aged expression and 
the aging depression A a , which indeed has the same func¬ 
tional form as for the subdiffusive CTRW and the hetero¬ 
geneous diffusion process. In the limit of strong aging, we 
also showed that ergodicity is seemingly restored and the 
disparity between ensemble and time averages becomes 
increasingly marginal. 

Confined aging SBM, in contrast, is qualitatively a 
quite unique process. Due to the time dependence of 
the diffusivity there is no thermal plateau for the 

ensemble averaged MSD. Instead this quantity is contin¬ 


uously decaying (subdiffusion) or increasing (superdiffu¬ 
sion). As shown here the functional behavior of confined 
aging SBM is remarkably rich. Concurrently, the time 
averaged MSD exhibits an intermediate plateau. In the 
presence of aging, we observe a deviation from the non¬ 
aged power-law behavior of the ensemble averaged MSD 
at longer times and a universal convergence to a plateau 
value. For strong aging, we again observe the conver¬ 
gence of ensemble and time averaged MSDs. We note 
that the behavior of confined CTRW is opposite: the 
time averaged MSD exhibits a power-law growth with 
exponent 1 — a, while the ensemble averaged MSD con¬ 
verges to the thermal plateau value [33]. 

In addition to these quantities we considered the first 
passage time density in the semi-infinite domain. In con¬ 
trast to the non-aging fractional Brownian motion we 
found a crossover between two characteristic scaling laws 
depending on the competition between aging and process 
time, t a and t. 

When using SBM as a stochastic model cognizance 
should be taken of the fact that it is a highly non¬ 
stationary process. The time dependence of its diffusiv¬ 
ity corresponds to a time dependent temperature (noise 
strength), and is therefore physically meaningless as de¬ 
scription for a system coupled to a thermostat. There 
exist, however, cases, in which SBM may turn out to be 
a physically meaningful approach. For instance, it was 
demonstrated that SBM provides a useful mean field de¬ 
scription for the motion of a tagged particle in a granular 
gas with a sub-unity restitution coefficient in the homo¬ 
geneous cooling phase & 

A number of aging features are quite similar between 
subdiffusive CTRWs ;3f|, heterogeneous diffusion pro¬ 
cesses [13], and aging SBM, as shown here. While frac¬ 
tional Brownian motion is ergodic, transient deviations 
from ergodicity and transient aging occur under confine¬ 
ment [3, [5(| • To reliably distinguish these models from 
another, it is therefore imperative to employ other diag¬ 
nostic stochastic quantities with characteristic behaviors 
for the respective processes [13, HI, HU ■ 
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